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Abstract 

The aim of this work is to study how the thermodynamic temperature is 
related to the known thermometers for nuclei especially in view of studying the 
microcanonical phase transition. We find within the MMMC-model that the 
" S-shape" of the caloric equation of state e* (T) which is the signal of a phase 
transition in a system with conserved energy, can be seen in the experimentally 
accessible slope temperatures T s i ope for different particle types and also in the 
isotopic temperatures Tue~Li- The isotopic temperatures Tn-He are weaker 
correlated to the shape of the thermodynamic temperature and therefore are 
less favorable to study the signal of a microcanonical phase transition. We 
also show that the signal is very sensitive to variations in mass of the source. 
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In this work we are interested in testing the different experimentally accessible thermome- 
ters for nuclei in order to understand which quantity is best related to their thermodynamic 
temperature. It is our purpose to show that it should be in principle possible to measure 
the specific signal of a microcanonical phase transition in an accurate experiment. 

The concept of phase transitions is usually discussed in connection with macroscopic 
systems. In such systems phase transitions are recognized by divergences in quantities 
like specific heat c(e*) = de* / dT t hd, where e* is the specific excitation energy and T t hd the 
thermodynamic temperature. If we are interested by similar phenomena in finite systems, the 
divergences are unsuitable to recognize and to classify phase transitions, since no divergences 
can occur. In a finite system the conservation laws become significant for the appearance 
and the shape of a phase transition. If a finite system is in thermal contact with a heat 
bath of constant temperature T t hd, i-e. is a canonical system, then the signal of a first order 
phase transition will become smeared showing up as a bump (instead of a divergence) in 
heat capacity or, equivalent, as a smooth anomaly in the caloric curve e*(T t hd). For a finite 
and isolated, i.e. microcanonical system the signal changes qualitatively. Here the total 
energy E of the system is a strictly conserved quantity. The first order phase transition 
becomes signaled by an "S-shape" in the caloric equation of state T t hd(E) as shown in 
the following figures in this publication. 

Since excited nuclei are an example of finite as well as isolated systems, they are especially 
suitable to study the signal of a microcanonical phase transition. This signal with an "S- 
shape" was frequently obtained in calculations, though it is still a subject of controversial 
discussion. Therefore it is of fundamental interest to test the theoretical findings by an 
experiment. 

Another important point is the fact that the thermodynamic temperature which can be 
simply measured in macroscopic physics is not accessible directly in a finite system. In a 
macroscopic system the measurement is simple through the fact that the size of the ther- 
mometer which gets into thermal contact with the system of interest is negligible compared 
to the size of the system. For a finite system this is obviously not satisfied. Especially for a 
microcanonical system it is not possible to bring a thermometer into thermal contact with 
a system without violating the strict energy conservation. Thus to obtain information on 
the thermodynamic temperature one needs to find an experimental observable which is not 
a temperature, but keeps information on the behavior of Tthd- 

Several suggestions for such observables were made. The candidates for "nuclear ther- 
mometers" are the slope temperatures from Maxwellian fits of energy distributions and the 
temperatures deduced from the isotopic ratios. In this work we are testing the quality of 
these "nuclear thermometers" as concerning their ability to reproduce the shape of the mi- 
crocanonical caloric equation of state {CES) T t hd(E). Another approach can be found in 
U . 

We concentrate on the microcanonical phase transition from evaporation to asymmetric 
fission which was predicted in the fragmentation of hot nuclei within the Berlin statistical 
fragmentation model MMMC ||. Similar signals were also predicted by other statistical 
models for nuclear fragmentation 0|| and also for atomic clusters @-@|. An experimental 
observation in ref. [0] seems also to support the existence of this phase transition. 

To study the "nuclear thermometers" we first obtain the signal of a phase transition 
in the T t hd{E) curve within the MMMC- model. Then we calculate the signals of the two 
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thermometers in question, the caloric curves of T s i ope (E) and T isotop i C (E) . 

Let us first briefly discuss the basics of the model. The strictly microcanonical MMMC- 
model H assumes a hot compound nucleus to be formed in a nuclear collision. After frag- 
mentation of this compound system the fragments remain coupled and exchange nuclei as 
long as they are in close contact. Consequently, the system is assumed to equilibrate statisti- 
cally shortly after the break-up. The volume which is accessed by the equilibrated fragment 
configuration is called the freeze-out volume. This means in terms of thermodynamics that 
the collection of all possible fragment configurations represents the maximum accessible 
phase space Q(E) for a given freeze-out volume. Q(E) is restricted by the valid conservation 
laws, which are the conservation of mass, charge, linear and angular momentum and of the 
total energy of the system, and also by the geometrical constraints. In the simulation the 
most important geometrical constraint is the size of the freeze-out volume, which is taken 
to be spherical. The radius rf of this volume, which is the only simulation parameter of the 
MMMC-model, is for the energy region of 1 to 4 MeV per nucleon at about r/ = 2.2A 1 / 3 fm, 
which corresponds to approximately 6 times the normal nuclear volume. When the frag- 
ments (which can be in excited states) leave this volume they may de-excite as they trace 
out Coulomb-trajectories. 

The output of the MMMC calculation is a collection of freeze-out configurations which 
are supposed to be representative for the entire phase space Q(E). The thermodynamic 
temperature T thd of these configurations is calculated by 

1 dS ds 



thd 



dE BE* 



where E is the total energy of the system, S = ks InQ(E) is the entropy, fee the Boltzmann 
constant and s = S/A and E* = E/A with A the mass of the decaying nucleus. The caloric 
curve T thd (E*) is plotted in all figures of this paper as a solid line with circles. 

We test the "nuclear thermometers" in two steps. First we plot the caloric curves 
T s i op e(E*) and T isotopic (E*) obtained from the MMMC-events after performing the Coulomb 
trajectories. Next we subject the calculated events to the software filter of the INDRA setup 
[TT| and plot the filtered caloric curves. 

For the INDRA-filter we assume the source velocity as 8.1 cm/ns, which is a typical 
quasi-projectile velocity measured with INDRA in mid-peripheral Xe + Sn collisions at 



50 A. MeV bombarding energy ||11|| . After the filtering we select the complete events for 
which the detected total charge and the total momentum is greater than 80% of the initial 
charge and of the initial momentum of the source, respectively. These events are used for 
constructing the filtered caloric curves. 

We start with the slope temperatures T s i ope for protons, deuterons, tritons, 3 He and alpha 
particles. The calculated kinetic energy spectra were fitted with the surface-evaporating 
Maxwellian source formula |I2| for every particle type: 



da {E kin - B) (-^) 

oc — — 2 e T ^ , (2) 



dE kin T slope 

where E^n is the center of mass kinetic energy of the particles, B the Coulomb barrier and 
Tsiope, which is the slope of the raw spectra, which is the desired slope temperature. We 
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calculate the slope temperature for protons, deuterons, tritons 3 He, alpha and also light 
IMFs. 

Figures [I], 0, [3], £| and |5] show the comparison of T thd (E*) and T s i ope (E*) for unfiltered (left 
plot) and INDRA-filtered (right plot) events for p, d, t, 3 He and a. Our most important 
finding in all these curves is that the slope temperatures resemble the general shape of 
T t hd{E*) before and after the filtering. 

The caloric curves for Li, Be and B (not shown here) which we have calculated only 
without filtering, repeat also the general shape of the phase transition despite of big error 
bars in T slope . 

The unfiltered T s i ope (E*) for p, d, t systematically achieve values lower then the ther- 
modynamic temperatures, while the values for 3 He and a are close to those of T t hd- We 
think that the last finding is just accidental. The unfiltered T s i ope are close for all particles 
with Z = 1 and for those with Z = 2, and the shift between them of ~ 0.3 MeV is due to 
the higher Coulomb repulsion. Here one can see in a very simple way that the result of the 
MMMC calculation cannot be described by a Maxwellian source with a unique temperature 
T s i ope (E*). Still we see that using the Maxwellian fit just as a recipe we obtain a useful tool 
to extract a pseudo-temperature which is correlated with the thermodynamic temperature. 

Let us now proceed to the isotopic temperatures JTSJ. The basic assumption for the 



isotopic temperature formula is a thermal equilibrium between free nucleons and composite 
fragments within a certain interaction volume V at a temperature T. The formalism is of the 
grandcanonical ensemble and ignores the possible effects of mass and energy conservation. 
The THe-Li(E*) isotopic temperature is given by 

Tue-u = lQ/ln(2A8 ^ Ll/ % Li ) , (3) 

1 3 He/ l^He 

and the Tn-He(E*) temperature by 

T^ He = U.3/ln(1.6 ^ H/ ^ H ) , (4) 

where Y(E*) is the particle yield. 

Fig. P shows a comparison of unfiltered and filtered TH e -Li{E*) with the thermodynamic 
temperature T thd (E*). Again we see that the signal of the phase transition survived the 
procedure. 

Performing the same for the H-He isotopic ratios, figure [7], we find that already the 
unfiltered T H _ He (E*) isotopic temperature is less sensitive to the "S-shape" in the T thd (E*) 
curve, but it also shows some structure at the phase transition. 

Finally we would like to address the question, how sensitive is the signal of the mi- 
crocanonical phase transition to the mass and charge of the compound nucleus. Figure |8| 
shows the T t hd(E*) curve for several masses and charges of the source. In the left panel we 
show that decreasing the charge of the source from Z=54 for 122 Xe (dots) to Z=50 for 122 Sn 
(empty triangles) does not influence much the phase transition. On the opposite, increasing 
the mass by 10 nuclei for 132 Xe (full triangles) shifts the transition signal to higher excitation 
energies by approximately 0.5 MeV per nucleon. 

In the right panel we show two additional curves for 80 Se and 250 Cf, thus strongly reducing 
and strongly decreasing the mass. Here the evaporation-fission phase transitions for 80 Se 
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appears as a very weak signal at E* « 1.5 A.MeV and T t hd ~ 3.8 MeV. The stronger signal 
at T t hd ~ 4.6 MeV is due to a different fission-multifragmentation phase transition. 

The fragmentation behavior of 250 Cf shows no phase transition at all, even for higher 
or lower E* . This is connected to the intrinsic instability against fission, so that even 
at very low excitation energies the nucleus fissions symmetrically instead of entering first 
the evaporation mode and later the asymmetric fission mode like the lighter nuclei do. 
This shows very plastically that there is no unique liquid-gas phase transition in nuclear 
fragmentation, but many different transitions depending on mass, and may be on other 
characteristics of nuclei. 

The dependence of the phase transition on the mass is natural, since increasing the 
mass we automatically increase the phase space Q(E*). Thus if changing A with excitation 
energy can also produce unusual shapes of the T t hd(E*) curveQ without undergoing any 
phase transition. This discussion serves mainly to make the point that to measure the 
microcanonical phase transition it is essential to keep especially the mass of the source at 
constant over the whole energy range. 

Further, it is important to obtain energy bining less than 0.2 MeV per nucleon to detect 



the discussed transition. The investigation in ref. |L6| misses the phase transition by per- 
forming energy steps AE* « 1 A.MeV. The current INDRA setup can realize AE* « 0.4 
A.MeV for the energy region of 1 to 4 A.MeV, which is at the limit of desired accuracy. If 
an additional experimental setup could measure the mass of the projectile fragments one 
could obtain much better E* reconstructions and AE* below 0.2 A.MeV. 

Summarizing the above we studied how the shape of the caloric curve T t hd(E*) at a 
microcanonical phase transition is correlated to the shape of different nuclear thermometers 
T s io P e(E*) and T isotopic (E*) . The signal of a microcanonical phase transition is an "S-shape" 
in the thermodynamic temperature. The slope temperatures for proton, deuteron, triton, 
3 He and alpha and the isotopic temperature Tne-Li{E*) reproduce this "S-shape" at correct 
excitation energies for INDRA-filtered and unfiltered events. The absolute values of this 
curves vary and do not coincide with the thermodynamic temperature. This means that 
the decrease of temperature with rising excitation energy at a first order microcanonical 
phase transition can be measured experimentally. This is of fundamental importance for the 
systematic study of phase transitions in finite systems. 

O.S. is grateful to GANIL for the friendly atmosphere during her stays there. This work 
was supported by IN2P3. 



We suppose together with [15,1(| and in opposite to [17] that the curve shown in ref. [18] is just 
the effect of changing the mass of the source. 
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FIGURES 




FIG. 1. Caloric curves T^diE*) and T s [ ope (E*) for protons. Left panel: MMMC calculation, 
T s i ope (E*) unflltered. Right panel: MMMC calculation, r siope (S*) filtered according to the INDRA 
setup. 
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FIG. 4. Like figure § but for 3 He 
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FIG. 5. Like figure [l], but for alpha. 



Source : Z=50, A=122, r Q =2.2 fm 
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FIG. 6. Caloric curves T t hd(E*) and Tn e -Li{E*) from the MMMC calculation. Tue-Li is 
shown unfiltered and filtered according to the INDRA setup. 
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Source : Z=50, A=122, r Q =2.2 fm 
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FIG. 7. Like figure 0, but for H-He isotopic temperature Tn-He(E* 
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FIG. 8. Signals of phase transitions in the thermodynamic temperature T t hd(E*) from MMMC 
calculations for different masses and charges of compound nuclei. 
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